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Recently, a diffusion Monte Carlo algorithm was applied to the study of spin dependent inter- 
actions in condensed mattes-. Following some of the ideas presented therein, and applied to a 
Hamiltonian containing a Rashba-like interaction, a general variational Monte Carlo approach is 
here introduced that treats in an efficient and very accurate way the spin degrees of freedom in 
atoms when spin orbit effects are included in the Hamiltonian describing the electronic structure. 
We illustrate the algorithm on the evaluation of the spin-orbit splittings of isolated C, Tl, Pb, Bi, 
and Po atoms. In the case of the carbon atom, we investigate the differences between the inclusion 
of spin-orbit in its realistic and effective spherically symmetrized forms. The method exhibits a very 
good accuracy in describing the small energy splittings, opening the way for a systematic quantum 
Monte Carlo studies of spin-orbit effects in atomic systems. 



PACS numbers: 



I. INTRODUCTION 



In the last few decades, Quantum Monte Carlo (QMC) 
techniques were successfully applied to a large number of 
systems in different fields of physics, quantum chemistry 
and materials science^—. The high accuracy obtained 
and the possibility of treating a relatively large number 
of particles within an affordable computational cost have 
certainly been the key ingredients for such a success. A 
fundamental property of QMC methods resides also in 
their scalability, which, due to the introduction and re- 
cent proliferation of massively parallel systems, makes 
a good case for a future applicability to large quantum 
many-body systems, up to now mainly investigated by 
means of Density Functional methods. Several years ago, 
an extension of the diffusion Monte Carlo (DMC) algo- 
rithm, employing an efficient recasting of the Green's 
function, was applied to many nucleon systems with in- 
teractions depending on the spin and isospin degrees of 
freedom^. This seminal paper has opened the way to 
studies of broad classes of systems in which spin related 
effects play a fundamental role. Most successful appli- 
cations so far attain the field of nuclear physics^. How- 
ever, in the last few years the method was extended to 
many electron systems, and in particular to the study 
of the two-dimensional homogeneous electron gas^i and 
parabolic quantum dots£. These models are often used 
to describe quasi two-dimensional nanostructurcs built 
on semiconductor heterojunctions, where the confining 
potential shows an asymmetry giving rise to a transverse 
electric field interacting with the moving electrons. This 
coupling can be described by an effective Rashba spin- 
orbit Hamiltonian. The next natural extension of the 
method in condensed matter application is the study of 



spin-orbit (SO) effects in atoms, molecules and solids. 
Preliminary estimations show that QMC methods should 
be able to provide the necessary accuracy and affordabil- 
ity for bringing to a completely new level the theoretical 
investigation in this area. The extension is highly non 
trivial, mostly due to the technical issues that might in 
principle limit the accuracy of a QMC calculation in pres- 
ence of a non-central interaction. 

One of the main issues concerning the DMC method is 
the need for trial wave functions, i.e. accurate approxi- 
mations to the exact solution of the Schroedinger's equa- 
tion. The availability of a good trial wave function repre- 
sents a crucial challenge both at the computational and 
at a more fundamental level. The trial wave function not 
only affects the computational efficiency of the DMC cal- 
culation, but for Fermionic systems, determines the accu- 
racy of the final result as a consequence of the necessity of 
applying the so-called "fixed-node" approximation"^—. 
An additional problem which comes to the forefront in 
the case of spin-orbit interaction is the necessity to deal 
with operators which do not necessarily commute with 
the rest of the Hamiltonian. 

Both these problems are usually tackled via a prelim- 
inary variational Monte Carlo (VMC) approach. Re- 
garding the first challenge, VMC is an important, effec- 
tive way of dealing with optimization of trial wave func- 
tions including almost arbitrarily complex many-body 
correlation s 12 ' 13 . Concerning the second issue, the DMC 
expectations of operators (denoted as O) not commuting 
with the Hamiltonian is known to be biased by errors 
depending on the trial wave function accuracy as well, so 
that efficient optimization is again an important way to 
alleviate this problem. In DMC, such errors stem from 
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the way DMC expectation values are computed: 

< O > D MC= < * T J°y > . (1) 

Here </>t represents the trial wave function while ip is the 
DMC ground state of the system. If [O, H] = this in 
principle introduces no bias since tp is an eigenstate of 
H. If this is not the case, the estimate is affected by a 
bias depending on <pT, and it is usually corrected to the 
leading order by using the formula 



< ip\0\ip >~ 2 < 4>t\0\^ > - < 4>t\0\4>t >, 



(2) 



where integrals are normalized and the last term corre- 
sponds to the VMC estimate of the operator O over <f>T- 
From these arguments it is possible to grasp how an ef- 
ficient VMC algorithm capable of explicitly dealing with 
spin degrees of freedom is of significant importance, not 
only as a robust and reliable numerical method itself, 
but also as for the development of the more sophisticated 
QMC algorithms. 

In this paper we present the theory, the Monte Carlo 
algorithm, and the results assessing the possibility of de- 
veloping wave functions containing the necessary corre- 
lation at the antisymmetric level (i.e. in the part of the 
wave function which is usually described by a Slater de- 
terminant) by means of a variational procedure. The 
many body correlations that are usually introduced in 
order to describe the short range effects of the Coulomb 
interaction are neglected on purpose for the C atom in 
order to make the peculiar aspects of the spin-dependent 
algorithm clear. Two body correlations are instead intro- 
duced in calculations relative to the more relevant cases 
of heavier atoms. The paper is organized as follows. In 
section II the VMC algorithm and the structure of the 
spinorial wave functions used in the calculations are illus- 
trated. Section III presents applications to the evaluation 
of spin-orbit splittings in C, Tl, Pb, Bi and Po. Section 
IV is devoted to conclusions. 



II. METHOD 



where R here represents the space coordinates for the 
system. 

On the other hand, for Hamiltonians containing a spin- 
orbit interactions, the same simplification is not correct. 
The trial wave function will be, in general, made up with 
Slater determinants containing single particle spinors and 
all the spin variables should be in this case taken into ac- 
count explicitly. Spinors in general have nontrivial forms, 
and the spin state of the single particles could change as 
a function of the space coordinates. This means that in 
general one needs to perform a summation over a com- 
plete spin basis: 



J s 



< j> T \R,S>< R,S\Q\^t > 



(5) 



Clearly, one way or another, one has to sample the rele- 
vant degrees of freedom in order to evaluate the impact 
of the spin on the expectation values. One possible way 
is to sum over the spin variables, i.e. sample the 2 N dis- 
crete space. 

What is proposed here is an alternative approach in 
which integration takes the place of the discrete sum, 
and by exploiting the VMC algorithm this can be per- 
formed with good computational efficiency. 
The trial wave function which we employ in calculations 
is a single Slater determinant of single-particle spinors. 
Let us first consider the simple two-electron case in order 
to illustrate the method. The spinors we will use are 



ij'i 



V'2 
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(0) 



The spacial coordinates of the two electrons will be f[ 
and f* 2 while the spin coordinates can be parameterized 
as follows: 



si 



sin(ai)e l 



cos(ai) J S2 V cos ( a 2) 
The Slater determinant can be written as 



sin(a 2 )e i<52 



(7) 



Variational Monte Carlo is a very efficient algorithm 
for evaluating expectation values of observables over a 
given trial wave function through the computation of in- 
tegrals of the kind 



< 1p T \0\lp T > 



(3) 



When one restricts to Hamiltonians or operators with no 
explicit spin dependence, it is not necessary to perform 
the summations over the spin degrees of freedom. These 
become static variables, and can simply be treated as 
labels. As a consequence, once a spin state and spin 
labels have been specified, the equation above can be 
rewritten as: 



J dRrARHrjR) 2 ^ 



(4) 



$(fi, s\,V2, s 2 ) = dot 



< si, n H'i > <si,n\ip2> 

<S2,f 2 |V>i> <s 2 ,r 2 \ip2> 



(8) 

A complete spin basis for this system is ttit4->4-t)44-! 80 
that, given a second wave function $' of the same form 
as <&, the following equality holds: 

< $|$' >=< $| tt><tt |$'> + < $| t4><tl I*' > 
+ < $| 4t><4t 1$' > + < *| U><U I*' > • (9) 

The terms involved in this summation are: 



<ttl*>=^i(ri)V»i(ra) 
<U\^>=^l(fi)^(fi) 

<ii|$>=^(n)^ 2 )- 



-Vi(r-i)^(r5) 
-Vl(ri)VH^) 



(10) 
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Writing $ explicitly as 



$(ri,si,r 2 ,s 2 ) = 
(V>i(ri)e Wl sin(ai) + ^ 2 (n) cos(ai)) • 

(^?(f 2 )e^ 2 sin(a 2 ) + ^(r 2 ) cos(a 2 )) 
(V? (ri)e i51 Bin(oi) + Vl(n) cos(ai)) • 

(^i(r 2 )e l52 sin(a 2 ) + ^(f 2 ) cos(a 2 )) 



(11) 



it is possible to demonstrate how the summation © can 
be exactly rewritten as an integral over the spin param- 
eters ai j2 and 5\ t 2 in the following way: 



< $1$' >= Const 



2tt 



$ (fi , si , r*2 , s 2 )* (n. , si, r 2 , 
daida2ddidS2 dridr 2 



that the same relation holds for any iV. This property 
holds due to the relations (fT3|) . which actually ensure a 
correct selection of all the necessary terms. The basic 
idea is that up (down) states of any particle will only be 
matched to the corresponding up (down) states of the 
same particle, and this precisely corresponds to projec- 
tion over a fixed spin state. Since all combinations are 
taken into account by the presence of both up and down 
states for all particles, we have accomplished the summa- 
tion over the entire basis set. 

All interesting obscrvablcs O which can be estimated 
with VMC, when acting on a Slater determinant, require 
at most the use of a linear combination of some new 
Slater determinants. This can be expressed as 
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In order to prove the above statement we take into con- 
sideration the terms corresponding to < $| tt><tt 
1$' >, since all other terms can be obtained in the same 
way. First of all we notice that 

da.i \ daj sin(a,) sin(ay) = (5ij27r 2 
o Jo 

doti I da.j cos(ai) cos(aj) = 5ij2ir 
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doti / duj sin(aj) cos(aj) = 0. (13) 
o Jo 

This means that only those terms with pairs of sines and 
cosines of the same angles will give a non zero contribu- 
tion. Using these relations we can select only the inter- 
esting terms: 

d ai da 2 [ip\{ri)*e- lSl sin(ai)V> 2 (r 2 )*e" i52 sin(a 2 ) 

-V?(ri)*e" 151 sin(ai)^ 1 1 (f 2 )*e- i52 sin(a 2 )] 
•[^(fije* 1 sin(ai)^ 2 (r 2 )e l52 sin(a 2 ) 
-Vf(ri)e l51 sin(a 1 )V 1 1 (r 2 )e l52 sin(a 2 )] (14) 

which give 

Const ■ [tiifiTtf^y^in)^ (r 2 ) 
+V 2 (r 1 )*V 1 1 (r 2 )*Vf(r 1 )Vi 1 (r 2 ) 

-^faJVifa) VM)^)] , (15) 

and this is precisely what one would obtain from the 
products of the first row of Eq. (fTO)) . As mentioned 
above, all the other terms of ([9]) can be obtained in the 
same way, proving the equality Eq. (|12|) . Furthermore, 
one can easily see how the integration over the variables 
Si is not necessary since these variables appear only as 
a phase, which exactly cancels out in the non- vanishing 
terms. Besides this, though the proof of (fT2j) was explic- 
itly given for the particular N = 2 case, it can be shown 



< R,S\0\$ >Y,< R ^ S \^ > 



(16) 



and for linearity the expressions we found can easily be 
applied. At this point, equation §5§ can be rewritten as 



/ dRdatp^(R,a)tpT{R,a) 



Oif> T (R,a) 



/ dRdaip^,(R,a)ip T (R,ct) 



(17) 



where the spin coordinates of the system S have been 
rewritten in terms of the parameters set a using ([7]). 
As mentioned earlier, the proposed algorithm does not 
represent the only possibility for spin summation. A pos- 
sible alternative could be that of a VMC-sampled sum 
over the possible spin states and this would also repre- 
sent an efficient algorithm. 



III. APPLICATIONS 

A broad class of problems exists in which the method 
discussed above can be applied. Nevertheless, it must 
also be stressed that it is not always easy to find a good 
trial wave function when an interaction of the spin-orbit 
kind is included in the Hamiltonian. This is mainly due 
to the non-locality of the spin-orbit interactions. 
In order to test the applicability of the method to atomic 
systems in presence of spin-orbit interactions, we chose 
to test this method on isolated C, Tl, Pb, Bi and Po 
atoms. By investigating both a light atom like Carbon 
and heavier elements, diverse spectra showing very dif- 
ferent energy splittings are considered, corresponding to 
the two distinct limits of LS and jj coupling. 



A. Carbon atom 

The carbon atom has six electrons and can be consid- 
ered as a light element. Although the spin-orbit effects 
are known to be very small and of little relevance for most 
properties, they still induce observable splittings in the 
energy spectrum. Furthermore, the computation of these 
very small energy differences represents a good starting 
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test for the efficiency of the method. The all-electron 
Hamiltonian which was employed in the calculations con- 
cerning the carbon atom is given by 



^ V 2m 



so — )+J2t~ 



(18) 



%<3 



where Vg accounts for the spin-orbit interaction for the 
i — th electron. It should be noted that the spin-orbit 
potential is the only rclativistic effect contained in this 
Hamiltonian. The rest of the relativistic corrections were 
neglected since they do not contribute to the level split- 
tings if we assume the wave functions given below. The 
trial wave functions employed in the calculations were 
linear combinations of Slater determinants of single par- 
ticle Hartree-Fock orbitals^. Since the interest, regard- 
ing the C atom, was focused on providing a good test for 
the method using a well-known type of wave function, 
no optimization was performed, and the combinations 
of Slater determinants were fixed by the imposed spin- 
spatial symmetries in the LS coupling. 
The spin-orbit interaction in atoms comes from an ap- 
proximated form of the Dirac equation and can be writ- 
ten (for the i — th particle) as 



Am 2 c 2 



(19) 



where V in this case is the total potential felt by the 
considered particle. This potential contains both the at- 
tractive contribution from the nucleus and the sum of 
the repulsive interactions with all other electrons. What 
is commonly done is approximating V with an effective 
potential V e ff having spherical symmetry. Within this 
approximation, equation (|19[) can be rewritten as 



ldV 



eff 



dr 



Li ■ Siip. 



(20) 



In case the Coulomb repulsion among electrons is ne- 
glected, one is left with a sum of single particle Hamil- 
tonians and if the factor in Eq. (f20f multiplying Li ■ Si 
is substituted by a constant, the problem can be exactly 
solved in analytic form. The term Li ■ Si in fact commutes 
with the rest of the Hamiltonian, which is spherically 
symmetric. However, as soon as the electron-electron 
Coulomb interaction is taken into account the exact so- 
lution is unknown. As already mentioned, since the SO 
interaction in C is known to yield very small energy split- 
tings, it is reasonable to treat it within the LS coupling 
procedure. This is done by writing the Hamiltonian (|18[) 
as H = Hi + H 2 where H 2 = J2i Kso an( i combining 
eigenstates of Hi (which are eigenstates of both L and 
S) in order to obtain eigenstates of J. Calling the eigen- 
states of Hi as \jLSMlMs > (7 stands for all resid- 
ual quantum numbers), it may be shown by using the 
Wigncr-Eckart theorem that 



< jLSM l Ms\H 2 \jLSM' l M' s > = 
A < 1 LSM L M S \L ■ S\jLSM' L M' s > 



where A is a constant depending on (-fLS). The eigen- 
states of J can be obtained as linear combinations 
of the [jLSMlMs > states and can be denoted as 
YjLSJMj >. At this point it is easily shown that 



< ~fLSJMj\H 2 \jLSJMj >= 
-A [J(J + 1) - L(L + 1) - S(S + 1)] 



(22) 



Notice how the expectation value of a sum of single par- 
ticle terms H 2 is related to that of L • S (S and L are the 
total spin and angular momentum) by the constant A. 
This constant contains the effect of the effective potential 
V e f f and in general it cannot be determined analytically. 
However, when the SO interaction (|19[) is approximated 
by Chi ■ Si with C constant, A can be calculated exactly. 
A first test for our method consisted in verifying the re- 
lations above and comparing the numerical result for A 
with the analytically computed value. In table U numeri- 
cal and analytical results are compared for three different 
states, showing excellent agreement with the predictions 
of the Wigncr-Eckart theorem. 



L S 


J 


numerical 


analytical 


numerical 


analytical 








> < £ . U ■ Si > 


A 


A 


1 1 





-1.001(2) 


-1.0 


0.5005(1) 


0.5 


1 1 


1 


-0.5004(7) 


-0.5 


0.5004(7) 


0.5 


1 1 


2 


0.4997(5) 


0.5 


0.4997(5) 


0.5 



TABLE I: Analytical and numerical results for the mean value 
of the SO interaction (C=l) and effective coupling constant A 
over the three lowest energy states given by fixed L = 1,S — 1 
and J = 0, 1, 2. Results are reported in atomic units 

A second set of calculations was then performed, in- 
cluding this time a SO interaction of the form (|2U|) in or- 
der to obtain an estimate of the energy splitting induced 
by the effective spherical potential V e f / . The only contri- 
butions to the SO splitting in the C atom come from the 
two occupied 2p orbitals, therefore it was only necessary 
to take the effective potential V e ff for 2p electrons into 
account. Following Slater—, in a single particle picture 
one could write the effective potential for the particle i 
as 



Zt ff (r) = Z- [ r j2\M*)\ a dr 



(23) 



(21) 



VMC calculations were done for the energy splittings of 
the lowest energy states (L = 1,S = 1). The results 
(see table |H]) are in qualitatively good agreement with 
experimental data^. The relative energy differences be- 
tween (J = 0,J = 1) and (J = 1,J = 2) states fol- 
low the theoretical ratio (1/2), though they both are too 
large. The fact that splittings from experimental data 
do not show the theoretical 1/2 ratio makes the result 
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for the larger splitting worse. In order to check whether 
our VMC method would give a correct description also 
in case of a realistic SO interaction, further testing was 
done: in case V e ff is substituted by Z/r 2 and the radial 
parts of the 2p single particle wave functions are given 
by a single orbital of the Slater type 



S fi (r) = Nj^- 1 exp(- Zjl r) 



= (2z jl p' +1 / 2 /[(2n jl )\} 1 / 2 



(24) 



(riij = 1,2,3,.. and Zij are fitting parameters) one can 
obtain an exact analytical form for the SO splitting. Also 
in this case the numerical results were in excellent agree- 
ment with the theoretical predictions. 



L S 


J 


realistic Vso 


effective Vso 


1 1 





-2.9(2) -10" 5 


-3.3(2) -10- b 


1 1 


1 


-1.4(1) -KT 5 


-1.6(1) -1(T 5 


1 1 


2 


1.4(1) -lcr 5 


1.7(2) -10" 5 



TABLE II: Numerical results for the mean value of the re- 
alistic (|19fl and effective SO interaction (|20[) per electron for 
the three lowest energy states. Results are reported in atomic 
units 

In order to check the effects of the approximations con- 
tained in the effective spherical potential, calculations 
were done substituting (j2U)) in the Hamiltonian with the 
realistic SO interaction (fH))) . Once more, SO potential 
expectation values were computed for the three lowest 
energy levels. In this case the energy splittings are re- 
duced with respect to the previous case becoming closer 
to experimental data (see Table III). The relative spac- 
ings retain the theoretical ratio 1/2 predicted for spher- 
ical potentials as before, but the effective SO coupling 
constant is close to the experimental one, at least for the 
smaller splitting (J = 0, J = 1). It remains however 
about 30 percent too large in the ( J = 1, J — 2) case. 
In the first case the effective screening of the inner elec- 
trons appeared to be not effective enough and one could 
be led to the conclusion that Hartree-Fock (HF) orbitals 
do not provide an accurate charge distribution. However, 
this last calculation showed how, at least for the lowest 
splitting, the screening effects affecting the SO coupling 
are reasonably well reproduced in the C atom already at 
the HF level. 



J VMC (effective Vso) VMC (realistic Vso) experimental 



-10(1.5) ■ 10" 

1 0.0 

2 20(1) ■ 10~ 5 



-9(1.5) • 10" 

0.0 
17(1) ■ HT e 



-7.5 ■ 10" 

0.0 
12.3 • 10 _ 



TABLE III: All electrons SO splittings (defined as E(J) - 
E(l)) for C. In the second and third columns VMC predic- 
tions relative to effective SO and realistic SO respectively. 
Experimental values in the fourth columns. Results are re- 
ported in atomic units 



B. Tl, Pb, Bi and Po atoms 

In order to give a more complete picture of the ap- 
plicability of the algorithm, the VMC method was also 
applied to the Tl (Z=81), Pb (Z=82), Bi (Z=83) and Po 
(Z=84) atoms, which exhibit sizable SO interaction ef- 
fects. While for C the LS coupling was employed, due 
to the increased strength of the SO interaction, a better 
approximation for the description of the energy splittings 
in heavier atoms is given by the jj coupling scheme. For 
this reason, the wave function is obtained by combining 
single-particle Ji (for the i — th particle) states into the 
cigenstates of the total angular momentum J. 
Due to the large number of electrons in the Tl to Po 
atoms, and the very deep ionic potential, it is customary 
to limit the explicit degrees of freedom to valence elec- 
trons by introducing an effective description in terms of 
pseudopotentials (or effective core potentials). Although 
there is no a-priori limitation in the number of electron 
that can be managed in a VMC calculation, very large 
energy fluctuations coming from the core states are very 
expensive to average out and are of marginal interest 
in determining the quantities we are interested in, given 
that SO splitting is only related to the presence of occu- 
pied valence p-states. The pseudopotential employed was 
taken from Kuechle et.alJ^, leading to a Hamiltonian (in 
atomic units) of the following form 



(1/2) 5>? +V m 



X/ r .. > 



(25) 



'</ 



where the indices i,j run over the valence electrons and 
V pp is the sum of a spin-orbit averaged ab initio pseu- 
dopotential V av , an SO operator V so and a term repre- 
senting effective charge of the pseudo-nucleus. 
The use of this pseudopotential has the additional bene- 
fit of making possible the comparison to the HF results 
reported in the same reference. In our calculations, the 
radial parts of s and p orbitals were represented by an 
expansion in four uncontracted gaussians. The single- 
particle orbitals optimization was performed with a vari- 
ational approach using an adapted implementation of the 
FIRE algorithmic. This simple algorithm is based on a 
molecular dynamics approach to the search in the vari- 
ational parameters space. The parameter evolution is 
driven by the forces induced by the total energy partial 
derivatives while appropriate cooling is introduced as pa- 
rameters motion is stopped whenever the velocities direc- 
tions are opposite to the forces. The forces, in turn, are 
computed during the VMC runs through a sampling of 
the wave function derivatives with respect to the param- 
eters. Given the form of the trial wave function ipx, the 
corresponding differentiation of the energy expectation 
according to a parameter a is given by 



< ^ T \H\^ T > 
a i,T ~ a <Vt|^t> 



(26) 
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and can be evaluated following the method also employed 
by Sorellai^, as 

< ip T \H\ip T >< ^T\d a \^T > , < 1pT\Hd a \lp T > 



our energy differences agree with the other two data sets 
essentially within the statistical error bars. 



(< vtI^t >y 



(27) 

Given the physical relevance of SO splittings in the heavy 
atoms here considered, the influence of correlation effects 
was investigated by introducing in the wave function a 
Jastrow factor— of the Padc' form: 



N 



J(R) = exp £V(r«) 



(28) 



being a a variational parameter. This factor introduces 
two body correlations, which in VMC calculations are 
usually necessary to obtain a realistic description of 
the correlation energy. In principle, the Jastrow factor 
should depend on the relative spin state of the electron 
pair. In presence of SO interactions inducing a spin ro- 
tation, this requirement would lead to a correlated wave 
function of the form Yii j f(ja{rij)<Ji ■ f?j, which would in 
turn imply the necessity of computing a sum over all the 
possible two electron spin states. This sum grows as 2^, 
and becomes very quickly unmanageable. For this rea- 
son we prefer, for the moment, to completely neglect the 
spin dependence in the two-body wave function, focusing 
on the gross effect of introducing the short range corre- 
lations induced by the Coulomb potential. We stress, 
however, that the introduction of correlation effects al- 
ready implies a huge improvement with respect to the HF 
method, keeping the computational cost at a reasonable 
level. 

In computing the SO splittings, the lowest states corre- 
sponding to the total angular momentum J eigenvalues, 
obtained from the multiple occupation of single particle 
p orbitals were considered for each atom. Comparing for 
instance Pb and C, despite the equal number of valence 
electrons, due to the diverse couplings (Jj and LS, re- 
spectively), these show different splitting patterns. In 
particular, for Pb the ( J = 0, J = 1) splitting is larger 
than the difference between (J = f , J = 2). 
VMC results are reported in table HVl together with the 
HF results of Kuechleil and the corresponding exper- 
imental values^. Calculations were also performed in 
absence of correlation, by removing the jastrow factor. 
Since in this case the trial function is given by the Slater 
determinant of spinors we can compare our results with 
the corresponding HF calculation. A good compatibil- 
ity is found. Interestingly, correlation effects on the SO 
splittings appear to be smaller than the statistical er- 
ror, suggesting that SO effects are mainly determined by 
the single particle properties of the wave function. From 
the comparison among the different data sets we confirm 
that our VMC calculation shows the correct ordering of 
the SO splittings among the states considered. Indeed, 





j 


VMC uncorrelated VMC+jastrow 


HF 


pynpTi m pyi t a 1 


Tl 














1/2 


0.0 


0.0 


0.0 


0.0 




3/2 


0.030(5) 


0.031(5) 


0.033 


0.035 


I u 














o 


0.0 


0.0 


0.0 


0.0 




1 


0.032(4) 


0.033(4) 


0.029 


0.035 




2 


0.044(5) 


0.045(5) 


0.046 


0.048 


Bi 














3/2 


0.0 


0.0 


0.0 


0.0 




3/2 


0.054(4) 


0.048(6) 


0.057 


0.052 




5/2 


0.074(5) 


0.070(4) 


0.078 


0.070 


Po 














2 


0.0 


0.0 


0.0 


0.0 







0.049(8) 


0.043(8) 


0.042 


0.034 




1 


0.080(6) 


0.077(5) 


0.071 


0.076 



TABLE IV: SO splittings for Tl, Pb, Bi and Po. The second 
and third columns contain uncorrelated and correlated VMC 
results respectively. The fourth and fifth columns report HF 
results by Kuechle and the relative experimental values^. 
States associated with zero splitting values are taken as a 
reference. Results are reported in atomic units 



IV. CONCLUSIONS 

We introduced an extension of the VMC method capa- 
ble to realistically describe spin-orbit splittings in heavy 
atoms. Calculations were tested first in the light C atom, 
and then extended to a set of heavier open p-shell atoms 
(Ti to Po). We demonstrated that the algorithm pro- 
vides correct evaluation of the underlying integrals and, 
remarkably, leads to sufficient accuracy for resolving the 
small SO energy differences. 

For the case of the carbon atom the small SO splitting is 
already reasonably well captured within a non-rclativistic 
Hartree-Fock wave function, although an even more ac- 
curate description would require the use of a better vari- 
ational wave function. Interestingly, the investigation of 
the effects of the spherical effective SO term showed an 
appreciable improvement when the more realistic version 
of the SO coupling operator was employed. 
Calculations of SO energy splittings carried on in heavier 
atoms required a wave functions based on jj rather than 
LS coupling. In these systems the possibility of a direct 
optimization of wave functions was investigated. The ef- 
fects of electron-electron correlations have been studied 
by introducing in the wave function an explicit, though 
simplified, two-body Jastrow factor. At this level dynam- 
ical quantum correlations seem not to have a huge effect 
on the SO splittings estimation. 

We believe that the obtained results are encouraging and 
the method can be applied to more complex systems 
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while retaining the efficiency and robustness of the VMC 
implementation. 
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